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INTRODUCTION  AND  THEORY 


I.  Purpose  of  the  Investigation 

The  solubilities  of  gases  in  liquids,  including  the  solubilities  of  noble 
gases  in  various  organic  solvents,  have  been  studied  by  many  investigators 
(1-7).    However,  investigations  of  the  solubility  of  gases  in  molten  salt  systems 
have  been  studied  only  recently.    These  investigations  were  begun  as  a  con- 
sequence of  the  production  of  xenon  and  krypton  as  by-products  in  high- 
temperature  nuclear  reactors.    The  presence  of  these  gases  in  the  reactor 
could  possibly  lead  to  undesirable  perturbations  in  nuclear  reactivity.  It 
was,  therefore,  desirable  to  know  the  solubilities  and  their  temperature 
dependences  for  such  gases  in  the  molten  fluoride  mixtures  used  in  the  oper- 
ation of  these  reactors. 

Fused  salt  systems  offer  a  radically  different  solvent  medium  for  the 
study  of  gas  solubilities.    Recent  studies  of  the  solubilities  of  noble 
gases  in  molten  fluoride  eutectic  mixtures  (8-10)  indicate  that  the  solu- 
bility of  a  gas  in  molten  salts  depends  to  a  large  extent  on  the  size  and 
polarizability  of  the  gas  molecule,  and  on  the  surface  tension  of  the  melt. 

The  purpose  of  this  work  is  to  present  new  data  for  gas  solubilities 
in  molten  salts.    The  work  includes:     (a)  solubility  measurements  of  nitrogen 
in  molten  sodium  nitrate  from  623°K  to  727°K  and  at  pressures  ranging  from 
lliO  atm.  to  U26  atm. ,  (b)  the  calculations  of  the  heat  and  entropy  of 
solution,  and  (c)  an  explanation  of  the  results  in  terms  of  the  size  and 
polarizability  of  the  gas  molecule,  and  the  surface  tension  of  the  melt. 

II.  Theoretical  Considerations 

The  solubility  of  a  gas  in  a  liquid  depends  on  the  nature  of  the  gas 
and  the  liquid,  and  on  the  pressure  and  temperature.    The  solubility  is 
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relatively  small  in  roost  cases,  and  is  consequently  governed  by  the 
thermodynamics  of  dilute  solutions  (11,  12).    One  of  the  general 
features  of  these  principles  is  that,  at  constant  temperature,  there 
is  a  linear  relation  between  the  equilibrium  pressure  of  the  gag  above 
the  solution  and  the  concentration  of  gas  dissolved  in  the  solution. 
This  linear  relation  between  concentration  and  pressure  is  known  as 
Henry's  law,  and  is  given  by 

m  -  K  P  1 
-     -p-  l- 

where  m  is  the  mass  of  gas  dissolved  in  a  unit  volume  of  solvent  in  equil- 
ibrium with  the  gas  phase  at  pressure  P  (at  constant  temperature),  and  Kp 
is  Henry's  law  constant  which  is  independent  of  the  pressure.    This  relation 
is  expected  from  kinetic  considerations,  since  compression  of  a  gas  tends  to 
increase  its  solubility.    Henry's  law  constant  depends  on  the  temperature 
and  the  nature  of  the  gas  and  the  solvent,  and  is  the  slope  of  the  solute 
concentration  versus  pressure  curve,  as  the  solute  concentration  approaches 
zero. 

Henry's  law  can  also  be  expressed  in  another  familiar  form.    Since  the 
mass  of  gas  dissolved  per  unit  volume,  m,  is  the  concentration  in  g.  cm. 
then  it  is  proportional  to  the  concentration  expressed  in  moles  of  solute 
per  liter  of  solution,  since  the  volume  change  on  solution  is  considered  to 
be  negligibly  small.    Also,  the  pressure,  P,  of  the  gas  in  equilibrium  with 
the  solution  is  proportional  to  the  number  of  molecules  in  the  gas  phase  in 
terms  of  moles  per  liter.    Therefore,  if  £d  represents  the  concentration  of 
the  dissolved  gas  in  moles  per  liter  of  solvent,  and  Cg  represents  the 
concentration  of  the  gas  phase  in  moles  per  liter,  then  Henry's  law  is  given 

by 
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where  Kc  is  usually  known  as  the  distribution  ratio  or  Ostwald's  coefficient 
of  solubility  (13) .    The  significance  of  this  distribution  coefficient  will 
be  discussed  in  later  sections. 

The  molecular  basis  of  Henry's  law  can  be  understood  from  the  fact  that 
in  dilute  solutions  the  solute  molecules  will  be  isolated  from  one  another 
and  surrounded  by  solvent  molecules.    Henry's  law  applies  as  long  as  the 
solute  molecules  never  approach  one  another  closely  enough  to  interact  sign- 
ificantly, because  the  escaping  tendency  is  then  proportional  to  the  number 
of  solute  molecules  in  the  fixed  amount  of  solvent.    The  intermolecular 
interaction  potential  between  non-polar  molecules,  or  between  rotating 
dipolar  molecules,  decreases  as  r    ,  where  r  is  the  intermolecular  distance. 
Thus,  the  interaction  energy  decreases  very  rapidly  with  increased  separation. 
There  are  two  principal  exceptions  to  these  conclusions:    (a)  if  the  solute 
molecules  are  very  large,  such  as  high  polymer  molecules,  then  the  probability 
that  the  solute  molecules  are  close  together  is  much  greater,  and  (b)  if 
the  solute  molecules  possess  electrical  charges,  the  intermolecular 
potential  decreases  only  as  r~l  or  less  (12). 

The  temperature  dependence  of  Henry's  law  constant  is  expressed  by  the 
Van't  Hoff  equation  (12) 

d(ln  K  )/dT  =  AHsoln./ET2  3. 
where  AHsoln.  is  the  heat  of  solution.    Integration  of  this  equation  yields 
log  K    -  -AHsoin./2.303  RT  +  constant  Lu 
Therefore,  a  plot  of  log  Kp  versus  l/T  should  yield  a  straight  line  whose 
slope  is  given  by  -AHsoim/2.303  B  ,  from  which  the  heat  of  solution  can  be 
obtained.    If  the  solubility  increases  with  temperature,  this  slope  will  be 
negative,  and  the  heat  of  solution  will  be  positive.    The  reverse  is  true 
when  the  solubility  decreases  with  temperature. 


The  interpretation  of  gas  solubility  by  Uhlig  (13)  expresses  the 
solubility  as  a  function  of  the  molecular  radius  of  the  gas  and  the  surface 
tension  of  the  solvent.    The  gas  molecule  is  assumed  to  be  spherical,  and 
the  surface  tension  of  the  solvent  is  assumed  to  be  the  same  for  molecular 
dimensions  as  it  is  for  the  macroscopic  bulk.    It  is  proposed  that  when  a 
gas  molecule  of  radius  r  enters  the  solvent,  it  produces  a  spherical- 
cavity  of  essentially  the  same  radius.    A  certain  amount  of  work  is  done  in 
producing  this  cavity,  since  any  increase  in  surface  area  of  the  liquid  is 
associated  with  a  definite  energy  change.    This  energy,  E-p  is  given  by  the 
increase  in  area,  due  to  creation  of  the  cavity,  multiplied  by  the  surface 
tension  of  the  solvent, f  , 

Ex  -  Uitr2Y  5. 
where  r  is  the  gas  molecule's  radius  in  Angstroms,  and Y is  the  surface 
tension  of  the  solvent  in  dynes  cm.~'L.    In  addition  to  this  latter  term, 
there  are  energy  terms  arising  from  solute-solvent  interactions,  which 
account  for  repulsive  and  attractive  forces  at  the  intermolecular  distances 
characteristic  of  molecular  separations  in  liquids.    The  nature  and  magnitude 
of  these  forces  depend  on  the  specific  properties  of  the  solute  and  solvent 
molecules,  which,  in  general,  present  a  complex  picture  too  difficult  to 
resolve  using  present-day  theories  of  molecular  interactions.    However,  a 
qualitative  picture  is  possible,  since,  for  example,  interactions  would  be 
expected  to  be  large  and  positive  between  gas  and  solvent  molecules 
containing  polar,  groups,  with  or  without  a  resultant  moment,  and 
relatively  small  and  negative  for  an  elementary  gas  dissolved  in  a  non- 
polar  solvent.    All  the  energies  in  this  category  are  termed  "interaction 
energies"  and  are  designated  as  E2.    Then  the  energy  change,  AE,  of  the 
system,  in  transferring  a  molecule  from  the  solvent  phase  to  the  gas  phase, 
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is  given  by 

AE  ■  Unr2Y   -  E2  6„ 
where  the  interactions  between  the  gas  phase  solute  molecules  and  the  gas 
phase  solvent  molecules  are  neglected,  since  these  are  small  except  at  high 
solvent  vapor  pressures.    The  Maxwell-Boltzmann  distribution  theorem  gives 
the  distribution  ratio  of  the  solute  molecules  between  the  two  phases  as 

Kc  =  Cd/Gg  =  exp  (-AE/kT)  7. 
where       and  Cg  are  the  concentrations  of  the  solution  and  gas  phase 
respectively,  and  k  is  the  Boltzmann  constant.    The  following  relation  is 
obtained  by  combining  eq.  6.  and  7. 

In  Kc  =  -liftr2Y  AT  +  E2/kT  8. . 

From  the  above  equation,  one  would  expect  that  if  the  interaction  energy 
remained  constant,  the  solubility  would  increase  with  decreasing  gas  molecular 
radius  and  decreasing  surface  tension  of  the  solvent.    Also,  for  gases  in 
which  the  solute- solvent  interactions  are  expected  to  be  small  (for  example 
with  the  noble  gases),  for  a  variety  of  solvents,  the  solubility  would 
depend  primarily  on  the  first  term  of  the  equation.    For  this  type  of  system 
the  plot  of  In  Kc  versus  the  surface  tension  of  the  solvent  should  be  linear 
and  have  a  slope  proportional  to  the  square  of  the  molecular  radius  of  the 
gas.    The  intercept  of  the  plot  would  be  a  measure  of  the  interaction  energy. 

From  the  solubility  data  for  several  gases,  including  02,  N2 ,  H2,  and 
CHj^,  in  a  variety  of  solvents,  Uhlig  found  that  the  interaction  energy 
was  approximately  constant  for  any  one  gas  and  independent  of  the  nature  of 
the  solvent  (13).    Also,  if  li1Tr2Y  >    E2,  the  solubility  of  the  gas 
increased  with  temperature,  and  when  hTTr2Y  ^  E2,  the  solubility  decreased 
with  temperature.    The  gas  molecular  radii  were  calculated  from  the  experi- 
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mental  plots  and  were  found  to  be  smaller  than  those  calculated  from  van 
der  Waals  constants.    This  indicated  that  the  estimated  sizes  of  the 
spherical  cavities  were  perhaps  too  small,  especially  when  the  thermal 
motions  of  the  gas  molecules  are  considered. 

In  the  simplified  model  of  Blander,  Grimes,  Smith,  and  Watson  (9), 
the  theoretical  distribution  ratio,  Kc,  for  an  inert  gas  in  solution 
equilibrium  with  a  fused  salt,  is  derived  by  equating  the  free  energy  of 
solution  of  the  gas  to  the  free  energy  of  formation  of  holes  the  size  of  the 
gas  molecules  in  a  continuous  fluid  having  the  same  surface  tension  as  the 
solvent.    Although  the  model  is  naive  it  yields  an  interesting  correlation 
with  experiment.    If  a  molecule  does  not  interact  with  the  liquid  in  which 
it  is  dissolved,  the  solubility  of  the  gas  can  be  estimated  from  the  free 
energy  changes  for  the  following  steps: 

(a)  The  gas  is  expanded  from  the  concentration  Cg  in  the  gas  phase 
to  .the  concentration       of  the  solution.    The  free  energy  change  for 
this  step  is  given  by 

AC^  -  RT  In  (Cd/Cg)  9o 

(b)  The  gas  atoms  are  contracted  to  "point  particles"  and  are  poured 
into  the  liquid.    The  free  energy  change  for  this  step  is  given  as 

AG2. 

(c)  The  "point  particles"  are  then  expanded  back  to  their  original 
sizes.    The  free  energy  change  for  this  step  is  given  by 

^3  =  k'AYmic  -AG2  10. 
where  A  is  the  area  of  the  hole  created  by  the  particle  when  it  is  expanded 
back  to  its  original  size,  Ymic.  is  the  microscopic  "surface  tension"  of 
the  solvent,  and  k«  is  a  constant  necessary  to  obtain  the  correct  units  of 
energy.    Since  the  molecules  are  assumed  to  be  spherical,  k'A  =  18. 08  r2, 
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where  r  is  the  effective  molecular  radius  in  Angstroms.    The  summation  of 
the  free  energy  changes  for  the  over-all  process  is  zero,  if  the  system  is 
in  equilibrium,  and  the  derived  relation  is  thus 

KC,theor.  =  £d/£g  "  e*P  <-18-08  £2T/RT)  11- 
where  is  generally  taken  as  the  macroscopic  value,  R  is  the  gas  con- 

stant, and  T  is  the  absolute  temperature.    This  model  neglects  solvation 
effects  in  comparison  to  the  large  endothermic  work  done  against  the  high 
surface  tension  of  a  typical  ionic  liquid  to  accommodate  a  gas  molecule. 
This  approach  is  justified  reasonably  well  by  the  cited  authors1  (9) 
observations  of  endothermic  heats  of  solution  for  various  noble  gases  in 
molten  fluoride  mixtures.    As  the  solvation  effects,  such  as  ion-induced 
dipole,  become  more  pronounced,  the  associated  exoenergetics  would  be 
expected  to  contribute  to  the  over-all  heat  of  solution  to  an  ever- 
increasing  extent.    At  some  point,  therefore,  the  solvation  energy  as- 
sociated with  a  given  solute  may  exceed  the  end oenerge tics  of  molecular 
cavity  creation,  resulting  in  the  more  usually  observed  exothermic  heat  of 
solution.    This  effect  would  be  even  more  enhanced  if  solvents  of  lower 
surface  tension  were  employed.    Increasing  solvation  effects  naturally  bring 
about  a  greater  gas  solubility  at  a  given  temperature  and  pressure,  which  in 
turn  causes  an  increasing  trend  in  the  ratio  Kc,exp. Ac, theor.    Such  was  the 
case  observed  in  reference  (9),  as  solvation  effects  increased  with  in- 
creasing polarizability  of  the  noble  gas  solute  atom  employed.    The  magnitude 
and  variation  of  the  solubility  with  size  of  the  gas  molecule  were  correctly 
predicted,  with  the  solubility  increasing  with  decreasing  size  of  the  noble 
gas  molecule.    Also,  the  endothermic  heats  of  solvation  increased  with 
increasing  atomic  weight  of  the  noble  gas.    The  theoretical  aspects  of  fused 
salt- gas  solvations  and  their  energetics  have  been  amply  reviewed  and  discussed 
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by  Stillinger  (lk) .    In  the  later  work  of  Shaffer,  Grimes,  and  Watson  OS), 
on  the  solubility  of  HF  in  molten  fluorides,  the  solvation  effects  obviously 
were  dominant  over  the  endothermic  work  against  surface  tension,  resulting 
in  fairly  high  solubilities  and  exothermic  heats  of  solution. 
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EXPERIMENTAL 

I.  Chemicals 

Reagent  grade  NaNC^  was  obtained  from  the  Baker  and  Adamson  Company, 
The  samples  were  heated  to  over  the  melting  point  to  remove  moisture. 
Residues,  if  any,  were  removed  with  a  porcelain  spatula,  or  the  salt  was 
decanted.    The  sample  was  cooled  in  a  desiccator,  and  was  then  crushed  to 
a  powder  for  weighing  into  a  Pyrex  bomb  liner. 

Nitrogen  gas,  99.9$%  pure,  bone  dry  grade,  was  obtained  from  the 
National  Cylinder  Gas  Company. 

II .  Apparatus 

A.    Commercial  Equipment 

A  500-ml  capacity  Inconel  metal  bomb,  Type  A2U3HC5  of  the  Parr 
Instrument  Company,  was  fitted,  according  to  specifications  by  J.L.  Copeland, 
with  two  CONAX  No.  TG-20-A-2-L  thermocouple  glands  packed  with  "Lava1"-,  a 
natural  magnesium  silicate.    Maximum  safe  working  pressure  of  the  bomb  was 
from  7,500  p.s.i.g.  at  2$°C  to  5,100  p.s.i.g.  at  600°C,  as  specified  by  the 
Parr  Instrument  Company  engineers.    A  Super-Gauge,  made  by  the  U.S. 
Gauge  Division  of  Ametek,  was  connected  to  the  bomb.    This  gauge  measures 
pressures  to  8,000  p.s.i.g.,  and  is  calibrated  in  divisions  of  100  p.s.i.g. 
The  accuracy  of  this  gauge,  as  quoted  by  the  manufacturer,  is  +  kO  p.s.i.g. , 
or  slightly  less  than  +  3  atm.  •  Both  of  the  CONAX  glands  contained  chromel- 
alumel  thermocouples,  one  of  which  measured  the  salt  temperature  and  was 
protected  by  an  Inconel  metal  tube  plus  an  outer  tube  of  Pyrex  glass.  The 
other  thermocouple  was  protected  by  a  Pyrex  tube,  and  was  used  to  measure 
the  gas  temperature.    The  manufacturer's  precalibrated  accuracy  of  each  of 
the  thermocouples  was  +  0.5°C    The  bomb  was  sealed  by  a  heavy-duty  copper 
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gasket  (seated  in  grooves  in  the  head  and  cylinder)  with  eight  cap  screws 
in  a  hardened  steel  ring,  split  in  halves.    These  ring  sections  are  retained 
by  a  heavy  steel  band.    A  photograph  of  the  entire  bomb  assembly  is  given  in 
reference  (16). 

The  vacuum  pump  used  was  a  Cenco  HYVAC  No.  91105.    When  a  gas  pressure 
greater  than  that  in  the  gas  cylinder  was  required,  an  Aminco  air-operated, 
single-ended  diaphragm  type  compressor  was  used.    This  vertical  unit  con- 
sists of  a  six-inch  single  shaft  air  operator  connected  directly  to  a 
hydraulic  plunger  which  transmits  pressure  to  a  hydraulic  medium.  This 
in  turn  causes  a  stainless  steel  diaphragm  to  flex  and  to  compress  gas  into 
the  head.    The  unit  is  rated  for  a  maximum  delivery  pressure  of  10,000 
p.s.i.g.    This  type  of  compressor  is  advantageous  in  that  the  gas  to  be 
compressed  comes  in  contact  only  with  a  metal  diaphragm  and  not  with  oil 
or  water. 

The  chrome 1-alumel  thermocouple  leads  from  the  bomb  were  connected  to 
a  Rubicon  potentiometer  No.  2703  with  an  internal  standard  cell.    A  -10° 
to  110°C  thermometer,  placed  at  the  potentiometer  terminals,  was  used  to 
record  room  temperature. 

B.    Constructed  Equipment 

A  Pyrex  liner  was  used  as  the  salt  container  to  protect  the  bomb 
from  corrosion.    The  diameter  of  the  liner  was  6  cm.  o.d.,  and  its  length 
was  If?  cm.    It  was  sealed  flat  on  the  bottom. 

The  furnace  was  constructed  by  W.C.  Zybko  (16).    It  consisted  of  one- 
quarter  inch  thick  Transite  outer  walls  and  a  U-inch  i.d.  Norton  Company 
aluminum  oxide  core  which  was  fixed  vertically  in  the  furnace  housing.  The 
outside  of  the  core  was  wound  with  18-gauge  Chrome 1-A  wire.    Norton  Company 
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Alundum  No.  RA  1055  refractory  cement  was  placed  over  the  outside  of  the 
wire-wound  core.    Super  Stic-Tite  asbestos  cement  served  as  insulation  between 
the  core  and  walls  of  the  furnace.    Temperature  in  the  furnace  was  maintained 
with  a  220- volt  Powerstat  and  a  0-5  amp.  a.c.  ammeter. 

Ill .  Procedure 

The  experimental  procedure  is  adopted  from  that  of  W.C.  Zybko  (16) „ 
The  reagent  grade  NaNO-^  was  melted  in  a  porcelain  casserole  to  remove 
moisture  and  other  impurities.    The  salt  was  then  transferred  to  a  desiccator 
and  allowed  to  solidify,  after  which  it  was  crushed  to  a  fine  powder  with  a 
mortar  and  pestle. 

The  Pyrex  liner  was  cleaned  with  dichromate  cleaning  solution,  and 
was  rinsed  with  deionized  water.    The  liner  was  dried  in  an  oven  at  about 
130°C  and  was  then  allowed  to  cool  in  a  desiccator.    It  was  then  charged 
with  30O-U5O  g.  of  the  dry,  powdered  salt,  weighed  to  the  nearest  0.5  g. 

Once  the  liner  was  in  the  bomb,  the  head  was  placed  in  position  and 
the  assembly  was  sealed.    Removal  of  air  and  any  additional  water  vapor 
was  assured  by  pumping  out  the  bomb  for  about  ten  minutes.    The  bomb  was 
charged  with  N2,  at  room  temperature,  to  some  desired  pressure.    The  com- 
pressor was  used  if  this  pressure  was  greater  than  that  in  the  gas  cylinder. 
The  bomb  was  allowed  to  stand  at  room  temperature  for  about  one  hour  to 
come  to  a  temperature-pressure  equilibrium.    After  equilibrium  measurements 
of  pressure  and  temperature  were  obtained,  the  bomb  assembly  was  heated  to 
an  initial  fused  salt  temperature  of  k5k°C  (except  in  one  run  when  the 
corresponding  pressure  involved  would  have  exceeded  the  prescribed  safety 
limits  of  the  bomb),  and  the  equilibrium  gas  pressure  and  temperature  were 
noted.    The  temperature  was  then  lowered  to  a  salt  temperature  of  Ui5°C, 
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with  the  system  being  agitated  frequently  to  ensure  equilibrium  and  to  avoid 
supersaturation  of  the  melt,  and  the  new  equilibrium  gas  pressure  and 
temperature  were  observed.    This  stepwise  procedure  was  continued  for  salt 
temperatures  of  U23°C,  1*13-5%  396. 5°C,  382.5°C,  368. 5°C,  and  355»5°C. 
This  experiment  was  repeated  five  times  for  the  same  eight  temperatures, 
each  time  with  a  different  starting  pressure  of       at  room  temperature. 
These  runs  yielded  eight  solubility  isotherms,  each  consisting  of  four  or 
five  solubility-pressure  points.    The  maximum  pressure  for  any  experiment 
was  1|26  atm.,  and  the  minimum  pressure  for  any  run  was  lUO  atm.  Molar 
volumes  of  N2  at  the  various  gas  temperatures  and  pressures  were  obtained 
from  graphs  of  the  P-^-T  data  compilation  of  Din  (17). 

As  in  previous  work  (18,  19),  the  decrease  in  moles  of  the  gas  phase 
in  the  system,  between  a  given  high  temperature  and  room  temperature,  was 
attributed  to  gas  solubility  in  the  melt.    The  density  of  liquid  NaNO-j  at 
a  given  temperature  and  pressure,  d-t,p>  used  for  calculating  the  pure  salt's 
total  volume,  was  obtained  from  the  following  empirical  equation  for  the 
density  in  the  absence  of  a  gas  phase 

dtjp  -  dtjl  +  (3.9  X  10"5)  P  g.  cm.-3  12. 
where  d^     is  the  density  at  t°C  and  1  atm.,  and  P  is  the  piston  pressure 
in  atm.    This  relationship  is  based  on  the  piston  pressure-density  data  of 
Owens  (20).    The  density  dt^ ^  is  obtainable  from  the  empirical  relation  of 
Bloom,  Knaggs,  Molloy,  and  Welch  (21) 

£t,i  =  2-13h  -  (7-03  x  10~k)  i  g-  cm-"3  13. 

A  minor  error  persists  as  a  result  of  using  pure  salt  densities  as  estimates 
of  the  actual  dilute  solution  densities.    It  is  thought  that  this  error  is 
within  the  experimental  error  of  the  procedure.    The  maximum  error  in  a 
pressure  measurement  is  slightly  less  than  +  3  atm.,  as  quoted  by  the 
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gauge  manufacturer  (18). 

Average  Henry's  law  constants,  Kp,  where 

Kp  =  £d/P    moles  of  gas  cm."^  atm.'^  lh* 
(P  is  the  gas-saturation  pressure  in  atm.),  and  distribution  coefficients, 
K  ,  were  calculated  for  each  solubility-pressure  isotherm. 


1L 


RESULTS  AND  CALCULATIONS 

The  determinations  of  the  solubilities  of  N2  in  fused  NaNO-^  at  the  various 
temperatures  and  pressures  were  based  on  the  difference  between  the  number  of 
moles  of  gas  in  the  gas  phase  at  room  temperature,  with  the  salt  in  the  solid 
state,  and  the  number  of  moles  of  gas  in  the  gas  phase  at  the  higher  temp 
erature,  in  equilibrium  with  the  molten  salt.    This  difference  in  the  number 
of  moles  of  gas  present  was  calculated  using  compressibility  data  for  Ng, 
which  was  obtained  from  detailed  plots  of  the  P-V-T  data  compiled  by  Din  (1?) . 
The  volume  of  gas  in  the  bomb  was  calculated  by  subtracting  the  volume  of 
the  salt  from  the  effective  volume  of  the  bomb  with  the  empty  liner.  The 
solid  salt  volume  was  calculated  using  a  density  of  2.261  g.  cm."^  (23). 
The  molten  salt  volumes  at  the  various  temperatures  and  pressures  were  cal= 
culated  using  densities  obtained  from  eq.  12  and  13.    Thermal  expansions  of 
the  bomb,  Pyrex  liner,  and  thermocouples  were  ignored  since  their  contri- 
butions to  the  total  bomb  volume  are  quite  small. 

Table  I  summarizes  the  experimental  pressures,  solubilities,  Henry's 
law  constants,  and  distribution  coefficients.    In  addition,  theoretical 
values  of  the  distribution  coefficients,  calculated  from  eq.  11,  are  also 
included.    Surface  tension  values  used  in  eq.  11  were  calculated  from  the 
empirical  equation  of  Addison  and  Coldrey  (22)  as  applied  to  fused  NaNO^ 

Y  =  116.6  -  0.000 (t  -  308)  dynes  cm."1  15. 
where  t  is  the  temperature  in  °C.    An  approximate  close  packed  radius  of 
2.00  2  was  taken  for  the  N?  molecule,  regarded  as  a  sphere  (21;).    It  is 
interesting  to  note  that  a  sample  calculation  for  the  expected  value  of  r  for 
N2,  using  eq.  11  together  with  an  experimental  value  of  Kc  for,  say,  6i±lo5°K, 
yields  a  value  of  r  =  1.16  S.    This  naive  approach  gives  a  value  which  appears 
to  be  too  small,  especially  when  Blander,  et.  al.,  (9)  were  forced  to  accept 
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effectively  larger  radii  than  usual  to  obtain  complete  agreement  of  theory 
with  experiment.    The  errors  recorded  for  the  average  Henry's  law  constants  are 
the  probable  errors  for  a  single  observation. 

Figure  1  is  a  graph  of  the  eight  solubility-pressure  isotherms,  with 
temperatures  as  indicated  in  °K.    The  plots  are  drawn  with  slopes  equal  to  the 
tabulated  average  Kp  values.    Individual  points  are  not  shown  for  the  sake  of 
avoiding  confusion;  they  are  available  for  graphing  from  Table  I.  The 
asterisked  values  of  Kp  in  Table  I  were  statistically  discarded  according  to 
Chauvenet's  principle  (25).    Figure  2  depicts  the  relationship  between  the 
common  log  Kp,av.  and  l/T  in  (°K)-1.    The  least  squares  straight  line  relation- 
ship for  this  graph  is 

log  Kp  -  (597  +  19)  (1/T)  -  (6.6U  +  0.03)  16. 
where  the  errors  are  the  least  squares  probable  errors,  and  units  Of  K  are 
those  of  Table  I.    The  theoretical  thermodynamic  equation  for  such  a  plot  is 
(15) 

lnKp«  -  -  (  AH/R)  (1/T)    +  ASp/R  17  „ 

where  Kp1  is  in  units  of  mole  l."1  atm."1,  AH  is  the  heat  of  solution,  and 
AS   is  the  entropy  of  solution  in  terms  of  equal  solute  gas  pressures  inside 
and  outside  of  the  liquid.    Use  of  eq.  16  and  17  yields  the  experimental  values 

AH  =  -2.73  +  0.09  kcal.  mole"1  and   ASp  =  -16.6  +  0.1  e.u., 
where  the  errors  are  again  probable  errors. 

Calculations  of  the  experimental  distribution  ratios  Kc>exp  ,  are  given 
as  defined  by  eq.  2.    It  is  assumed  that       is  small  enough  to  render  the 
solution  ideal,  so  that  the  activity  of  the  dissolved  gas  may  be  equated  to 
its  concentration.    Furthermore,  ideal  gas  behavior  is  assumed  for  Ng  at  the 
high  temperatures,  in  spite  of  the  high  pressures  involved,  so  that  the 
fugacity  may  be  equated  to  the  pressure. 
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Figure  1.    Solubility-Pressure  Isotherms  for  Nitrogen  in  Fused  Sodium 
Nitrate.     (Temperatures  are  those  of  the  Fused  Salt  in  °K) . 
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Figure  2.    Temperature  Dependence  of  Henry's  Law  Constant  for  Nitrogen 
in  Fused  Sodium  Nitrate. 
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DISCUSSION 

The  negative  temperature  coefficient  of  solubility  for  N2  in  fused  NaNO^, 
and  the  attendant  exothermic  heat  of  solution,  lead  to  the  conclusion  that 
solvation  effects  in  this  system  are  of  greater  importance  than  the  work 
required  to  create  a  solute  molecule  cavity.    This  conclusion  is  also  supported 

by  the  large  lCfexpt/Kc>fheor  ratios>  Ic,expt(avera^e)  bein&  from  23  to  a3 
times  greater  than  KC)theor*    ^^-s  rati°  *s  a^so  seen  to  decrease  with  temper- 
ature.   However,  no  immediate  significance  can  be  attached  to  this  trend. 
Prior  work  on  solubilities  of  the  nobla  gases  in  molten  fluorides  has  usually 
resulted  in  this  ratio  being  less  than  unity,  with  the  smallest  values  for  the 
less  polarizable  gases  such  as  He  and  Ne  (9,  10).    In  some  of  these  cases,  with 
more  polarizable  gases  such  as  Ar  and  Xe,  the  ratio  was  sometimes  slightly 
greater  than  unity,  presumably  owing  to  the  effect  of  greater  ion-induced  dipole 
interactions  in  augmenting  the  solubility. 

Results  from  the  work  on  solubility  of  HF  in  molten  fluorides  gave 
experimental  heats  of  solution  ranging  from  -3.85  kcal.  mole-1,  for  h$  mole  % 
of  NaF  in  a  NaF-ZrF^  mixture,  to  -9.70  kcal.  mole"1  for  80.5  mole  %  NaF  (15) . 
The  observed  AR  of  -2.73  kcal.  rnole-1  for  N2  in  NaNO^  is  somewhat  smaller  than 
the  lowest  such  value  for  the  HF  work,  implying  less  exothermic  solvation 
effects  for  the  N2  system  than  for  the  HF  systems.    This  is  as  expected  in 
view  of  the  exoenergetic  contribution  of  ion-permanent  dipole  interactions 
which  must  occur  in  the  latter  systems.    Also,  the  Henry's  law  constants 
observed  in  the  HF  work  are  about  ten  times  greater  than  the  values  observed 
for  the  N2  system.    This  may  again  be  partially  attributable  to  such  ion- 
dipole  solvation  effects  occurring  in  the  HF  systems.    In  addition,  the 
relatively  high  stability  of  NaF-HF  complexes  occurring  in  these  systems  no 
doubt  serves  to  increase  HF  solubility  and  to  make  the  heat  of  solution  still 
more  negative  (1$) . 
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The  AS    entropies  of  solution  found  in  the  HF  work  varied  from  ~lh<0  e.u., 
for  the  US  mole  %  NaF  system,  to  -15.2  e.u.  for  the  80.5  mole  %  NaF  system 
(15).    The  value  of  -16.6  e.u.  obtained  for  N2  in  NaNO-j  therefore  appears 
reasonable  in  comparison  to  these  values  for  other  diatomic  molecules.  These 
entropies  of  solution  for  N2  and  HF  diatomic  molecules  are  noticeably  more 
negative  than  values  obtained  for  monatoraic  noble  gases  in  fluoride  melt; 
(8,  9).    This  observation  was  predicted,  since  at  higher  temperatures  a  non- 
spherical  gas  molecule  rotates  more  freely  and  the  "surface  area"  of  its 
cavity  is  correspondingly  increased,  leading  to  a  more  negative  entropy  of 
solution  than  for  a  spherically  symmetric  rare  gas  atom  (8). 

Copeland  and  Zybko  have  observed  Henry's  law  constants  for  high  pressure 
solubilities  of  Ar  (18)  and  He  (19)  in  molten  NaNCL  at  61i2°K.  For 
Ar: 

£p,Ar  *  0-1.2  +  1.7)  x  Kf7  mole  cm.-^atm."1 
and  for  He: 

KpjHe  =  (22.7  +  0.7)  x  10"  7  mole  cm.'^atm"1. 
The  close-packed  "spherical*1  radius  of  the  N2  molecule  is  very  nearly  the 
same  as  that  of  Ar  (1.92  8)(9).    Thus,  the  cavity  creation  work  done  against 
the  comparatively  small  surface  tension  of  NaNO^  by  either  gas  would  be 
expected  to  be  small  and  approximately  the  same  for  both.    However,  the 
polarizability  of  N2  is  1.73U  x  10"2^  cm.-3  molecule-1,  as  compared  to 
1.626  x  lO"2^  cm. 3  molecule-1  for  Ar  (26).    On  this  basis,  one  would  expect 
slightly  greater  ion-induced  dipole  solvation  effects  for  N2  than  for  Ar, 
which  would  tend  to  make  the  former  slightly  more  soluble  if  all  other  effects 
are  nearly  equal.    This  is  demonstrated  by  the  slightly  greater  solubility  of 
N2  at  61|1.5°K.    The  polarizability  of  He  is  only  about  0.2036  x  10~2^  cm.3 
molecule-1  (26),  and  its  close-packed  radius  is  about  1.22  &  (9). 
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Thus,  much  less  work  against  surface  tension  is  required  to  create  a  hole  for 
a  He  atom  than  for  either  N2  or  Ar.    However,  solvation  effects  for  He  would 
be  practically  nonexistent.    It  would  appear,  therefore,  that  the  slightly 
greater  solubility  of  He  may  be  primarily  the  result  of  the  relative  ease  of 
introducing  such  a  small  atom  into  the  liquid,  in  spite  of  the  partially 
compensating  solvation  effects  which  exist  for  the  larger  solute  molecule*,: 


2h 


ACKNOWLEDGMENTS 

The  author  wishes  to  express  his  sincere  appreciation  to  Dr.  James  L„ 
Copeland  for  his  guidance  and  encouragement.    Without  his  assistance  this 
work  and  the  resulting  paper  would  not  have  been  published.  Appreciation 
is  also  extended  to  Dr.  Walter  C.  Zybko  for  his  construction  of  the  furn  ce 
necessary  to  accomplish  this  work. 

The  author  gratefully  acknowledges  support  of  this  work  by  the 
National  Science  Foundation,  Grant  GP-U27U. 


2? 


LITERATURE  CITED 


1.  A.  E.  Markam  and  K.  A.  Kobe,  Chem.  Rev.,  £8,  519  (19Ul) . 

2.  J.  H.  Hildebrand  and  R.  L.  Scott,  "The  Solubility  of  Non-electrolytes,4' 
3rd  ed.,  Dover  Publications,  Inc.,  New  York,  N.Y.,  1961a,  [488pp. 

3.  M.  W.  Cook,  U.  S.  Atomic  Energy  Commission,  UCRL-2l*59  (1951) . 

U.    H.  L.  Clever,  R.  Battino,  J.  H.  Sayler,  and  P.  M.  Gross,  J.  Phys,  C^em., 
61,  1078  (1957). 

5.    H.  L.  Clever,  J.  H.  Sayler,  and  P.  M.  Gross,  ibid.,  62,  89  (1958). 


6.    L.  W.  Reeves  and  J.  H.  Hildebrand,  J.  Am.  Chem.  Soc,  79,  1313  (1957). 


7.  E.  Nussbaum  and  J.  B.  Hursh,  J.  Phys.  Chem.,  62,  81  (1958). 

8.  W.  R.  Grimes,  N.  V.  Smith,  and  G.  M.  Watson,  ibid.,  62,  862  (1958). 


9.    M.  Blander,  W.  R.  Grimes,  N.  V.  Smith,  and  G.  M.  Watson,  ibid.,  63 
1161*  (1959).  ~ 

10.    R.  B.  Evans  III,  W.  R.  Grimes,  N.  V.  Smith,  and  G.  M.  Watson, 
J.  Chem.  Eng.  Data,  7,  285  (1962). 


11.  G.  M.  Barrow,  "Physical  Chemistry,"  McGraw  Hill  Book  Co.,  Inc., 
New  York,  N.Y.  I96I,  pp.  1*81-1*83. 

12.  G.  N.  Lewis  and  M.  Randall,  "Thermodynamics,"  2nd  ed.,  revised  by 

K.  S.  Pitzer  and  L.  Brewer,  McGraw  Hill  Book  Co.,  Inc.,  New  York,  N»  Y„, 
1961,  pp.  231-21*0. 

13.  H.  H.  Uhlig,  J.  Phys.  Chem. ,  1*1,  1215  (1937). 

11*.    F.  H.  Stillinger,  Jr.,  in  "Molten  Salt  Chemistry,"  M.  Blander  editor, 
Wiley-Interscience,  New  York,  N.  Y.,  1961*,  pp.  87-105. 

15.  J.  H.  Shaffer,  W.  R.  Grimes,  and  G.  M.  Watson,  J.  Phys.  Chem. ,  63, 
1999  (1959).  T" 

16.  W.  C.  Zybko,  "Effects  of  Inert  Gas  Pressure  and  Solubility  on  Electrical 
Conductance  of  Fused  NaNO-,,"  Ph.D.  Thesis,  Kansas  State  University,  1965, 
p.  15.  5 

17.  F.  Din,  "Thermodynamic  Functions  of  Gases,"  vol.  3,  Butterworth  and  Co., 
Ltd.,  London,  1961,  pp.  11*6- 150. 

18.  J.  L.  Copeland  and  W.  C.  Zybko,  J.  Phys.  Chem. ,  6£,  3631  (1965). 

19.  J.  L.  Copeland  and  W.  C.  Zybko,  ibid.,  70,  181  (1966). 

20.  B.  B.  Owens,  private  communications  to  J.  L.  Copeland,  1965. 


26 


21.  H.  Bloom,  I.  W.  Knaggs,  J.  J.  Molloy,  and  D.  Welch,  Trans.  Faraday  Soc, 
h9.  Hi58  (1953). 

22.  C.  C.  Addison  and  J.  M.  Coldrey,  J.  Chem.  Soc,  1961,  U68. 

23.  ♦•Handbook  of  Chemistry  and  Physics,*1  Chemical  Rubber  Publishing  Co., 
Cleveland,  Ohio,  Uoth  ed.,  1959,  p.  656. 

2k.    W.  J.  Moore,  "Physical  Chemistry,*1  Prentice-Hall,  Inc.,  Englewood  Cliffs, 
N.J.,  3rd  ed.,  1962,  p.  229. 

25.  J.  W.  Mellor,  "Higher  Mathematics  for  Students  of  Chemistry  and  Physics," 
Dover  Publications,  New  York,  N.  T.,  19U6,  pp.  563,  623. 

26.  E.  A.  Moelwyn-Hughes,  "Physical  Chemistry,"  Pergamon  Press,  Oxford, 
2nd  ed.,  1961,  p.  383. 


27 


VITA 

The  author  was  born  in  Winnsboro,  South  Carolina  on  January  22,  19U1. 
At  the  age  of  five,  he  moved  to  Greensboro,  North  Carolina,  where,  after 
attending  Mt.  Zion  Elementary  School,  he  received  his  high  school  education 
at  the  James  B.  Dudley  High  School.    He  enrolled  at  the  Agricultural  and 
Technical  College  of  North  Carolina  in  the  fall  of  1959  and  received  the 
degree  Bachelor  of  Science  in  June,  1963. 

He  entered  the  Graduate  School  of  Kansas  State  University  in  September, 
1963,  and  in  June  of  1966  he  was  awarded  the  degree  of  Master  of  Science  by 
Kansas  State  University. 

The  author  is  a  member  of  Beta  Kappa  Chi  National  Scientific  Honor 
Society,  Alpha  Kappa  Mu  National  Honor  Society,  and  Alpha  Chi  Sigma  National 
Chemical  Fraternity. 

Publications 

1.  C.  H.  Edwards,  L.  Seibles,  et.al.,  "Chromatography  of  the  Methyl  Esters 
of  Selected  Fatty  Acids  on  ^Tass  Fiber  Paper,"  J.  Chromatog.  II,  3U9-35U 
(1963).  ' 

2.  J.  L.  Copeland  and  L.  Seibles,  "The,  High-Pressure  Solubility  of  Nitrogen 
in  Fused  Sodium  Nitratej  Temperature  and  Pressure  Dependences,  and  the 
Heat  and  Entropy  of  Solution,"  J.  Phys.  Chem. ,  70,  (accepted)  (1966)0 


SOLUBILITY  OF  NITROGEN 
IN  FUSED  SODIUM  NITRATE 

by 

LAWRENCE  SEIBLES 

B.S.  (Honors)  The  Agricultural  and  Technical  College 
of  North  Carolina,  1963 


AN  ABSTRACT  OF  A  THESIS 

submitted  in  partial  fulfillment  of  the 

requirements  for  the  degree 

MASTER  OF  SCIENCE 

Department  of  Chemistry 

KANSAS  STATE  UNIVERSITY 
Manhattan,  Kansas 

1966 


The  purpose  of  this  work  is  to  present  new  data  on  the  solubility  of  gases 
in  molten  salt  systems.    A  Parr  Inconel  metal  bomb  with  chrome l~alumel  thermo- 
couples was  used  to  study  the  temperature  and  high-pressure  dependences  of  the 
solubility  of       in  molten  NaNO^. 

Eight  solubility-pressure  isotherms  for  N2  in  fused  NaNO^  were  obtained, 
each  consisting  of  four  or  five  solubility-pressure  points.    The  solubilicies 
ranged  from  2.75  x  10" ^  moles  cm"3,  at  6££. £°K  and  118  atm. ,  to 
7.I4O  x  10"^  moles  cm."3  at  686. £°K  and  lj26  atm.    Average  Henry's  law  constants, 
Kp,  varied  from  (20.2  +  0.I4)  x  10"  7  mole  cm."3  atm.-1  at  628. £°K  to 
(15.2  +  1.2)  x  10"7  mole  cm.-3  atm."1  at  727. 0°K.    Experimental  values  of 
AH  =  -2.73  +  0.09  kcal.  mole"1  and  ASp  =  -16. 6  +  0.1  e.U.  were  obtained 
from  the  plot  of  log  Kp  versus  l/T  and  the  theoretical  thermodynamic  relation. 

The  negative  temperature  coefficient  of  solubility  for  N2  in  fused 
NaNO^,  and  the  attendant  exothermic  heat  of  solution  are  interpreted  on  the 
basis  that  the  solvation  effects  in  the  system  are  of  greater  importance  than 
the  work  required  to  create  a  cavity  for  a  solute  molecule.    This  interpretation 
is  also  supported  by  the  fact  that  experimental  distribution  ratios.  Kc,  are 
larger  than  the  theoretical  values.    The  results  also  indicate  that  for 
molecules  of  nearly  the  same  size,  greater  solubility  in  a  given  medium 
occurs  for  the  more  polarizable  molecule. 


